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Geochemistry  of  Freezing  Brines 

Low-Tempera  tore  Properties  of  Sodium  Chloride 

VALERIE  L.  THURMOND  AND  GARRETT  W.  BRASS 


INTRODUCTION 

Brines  that  form  by  evaporation  and  mineral 
dissolution  in  regions  that  are  cold  or  arid  or  both 
are  common  and  have  been  the  subject  of  many 
geochemical  investigations  because  of  their  un¬ 
usual  chemical  properties.  We  are  unaccustomed 
to  thinking  of  brine  formation  and  mineral  precip¬ 
itation  by  freezing,  although  this  is  an  important 
process  in  cold  regions.  Some  brines,  particularly 
those  with  a  high  concentration  of  calcium  chlor¬ 
ide,  can  remain  entirely  unfrozen  at  -50°C  such  as 
those  in  the  dry  valleys  of  Antarctica.  Saline  solu¬ 
tions  in  subsea  permafrost  are  derived  from  sea¬ 
water  and  concentrated  by  freeze-thaw  processes 
and  reactions  with  soil  minerals.  In  frozen  ground 
the  unfrozen  water  due  to  the  presence  of  salts 
changes  the  frost-heaving  characteristics  of  soil 
and  reduces  the  stability  of  structures  such  as 
buildings,  roads  and  drilling  platforms.  Brine  in¬ 
clusions  affect  the  properties  of  sea  ice  and  the 
weathering  of  rocks,  soils  and  building  materials 
in  cold  climates. 

These  considerations  make  the  phase  relation¬ 
ships  in  freezing  brines  important  to  oceanog¬ 
raphy  and  cold  regions  studies.  Other  than  several 
experimental  studies  of  freezing  seawater  (Nelson 
and  Thompson  1954,  Thompson  and  Nelson  1956, 
Richardson  1976)  and  the  binary  and  ternary  solu¬ 
bility  data  in  Linke  (1958),  little  is  known  about 
the  changes  in  phase  composition  that  occur  when 
brine  mixtures  are  frozen.  Calculating  these 
changes  is  difficult  since  the  thermodynamic  prop¬ 
erties  of  aqueous  electrolytes  below  0®C  are  virtu¬ 
ally  unknown,  with  the  exception  of  some  specific 
heat  data  for  calcium  and  magnesium  chloride 
(D’Ans  et  al.  1977).  Thermodynamic  quantities  of 
interest  to  solubility  studies,  e.g.  activity  coeffi¬ 
cients,  enthalpies  and  heat  capacities,  change  rap¬ 


idly  with  concentration,  especially  at  low  tempera¬ 
tures,  yet  there  are  very  few  measurements  below 
0°C. 

Recent  improvements  in  the  calculations  of 
thermodynamic  properties  of  mixed  electrolytes  at 
room  temperature  and  above  by  Pitzer  and  his  co¬ 
workers  have  beat  applied  to  the  calculation  of 
mineral  solubilities  in  concentrated  brines  at  25  °C 
(Eugster  et  al.  1980,  Harvie  and  Weare  1980). 
These  results  shed  light  on  old  problems  such  as 
the  mismatch  between  the  calculated  or  experi¬ 
mentally  determined  solubility  of  salts  and  the  ob¬ 
served  sequences  of  minerals  deposited  by  the 
evaporation  of  seawater.  The  model  is  now  being 
extended  to  high  temperatures  and  pressures  by 
various  workers  for  use  in  geothermal  studies,  pet¬ 
roleum  engineering  and  mining  geology  (particu¬ 
larly  hot-brine  ore  deposition).  Up  to  now,  how¬ 
ever,  temperatures  below  25  °C  have  not  been  in¬ 
cluded  in  these  studies. 

The  purpose  of  this  work  is  to  extend  the  treat¬ 
ment  of  electrolytes  to  freezing  temperatures  (-40 
°C),  beginning  with  the  properties  of  sodium 
chloride.  Sodium  chloride  is  a  major  constituent 
of  many  freezing  brines,  and  the  thermodynamic 
properties  of  its  solution  at  room  temperature  are 
well  known.  We  have  measured  heat  capacities  of 
sodium  chloride  solutions  from  25  ®C  to  -40  °C 
with  a  differential  scanning  calorimeter.  Combin¬ 
ing  these  results  with  literature  values  at  25  °C,  we 
have  calculated  the  activity  coefficients  from  0.5 
molal  to  saturation  and  down  to  -40  °C  using  Pit- 
zer’s  equations.  The  solubility  of  ice  (freezing 
point  curve)  is  accurately  predicted  from  the  low- 
temperature  activities  of  the  saturated  solutions. 
The  success  of  the  method  indicates  that  activities 
of  other  salts  can  be  determined  at  very  low  tem¬ 
peratures  and  that  these  can  be  used  to  predict  the 
properties  of  mixtures.  Work  is  underway  to  ex- 


lead  the  calculations  to  the  seawater  system  Na- 
Ca-Mg-K-Cl-SO*. 


METHODS 


The  solubility  product  (K°)  is  a  function  of  tem¬ 
perature  and  pressure  only  and  for  sodium  chlor¬ 
ide  is  given  by 

K %T,  1  atm)  -  «k,Cf'>'±N*a »tm)  (1) 
where  m  =  molality 

y±  »  activity  coefficient  at  1  atm  pressure 
T  «  temperature. 

If  the  temperature  variation  of  y±  is  known,  the 
equilibrium  constant  can  be  calculated  at  any  tem¬ 
perature  from  and  the  experimental  solubility 
(m).  The  temperature  dependence  of  y±  over  a 
wide  temperature  range  can  be  found  by  examin¬ 
ing  the  temperature  dependence  of  the  partial 
molal  free  energy  (a6°)  or  the  excess  free  energy 
(AG").  The  different  methods  have  been  reviewed 
by  MiDero  (1979). 

The  recent  work  of  Silvester  and  Pitzer  (1977) 
and  Rogers  and  Pitzer  (1961)  demonstrates  the  use 
of  excess  free  energies  for  determining  the  temper¬ 
ature  dependence  of  y±  and  the  activity  of  water 
(aw).  The  total  excess  free  energy  (C“)  is  given  by 

G**  mn*  »mRT( l-$+lny±)  (2) 

where  nw  -  mass  of  water  (kg) 

w  m  number  of  ions  per  formula  unit 
m  -  molality  (mol  kg*1) 

4>  •  osmotic  coefficient 
R  m  gas  constant 
T  -  temperature  (K). 

The  osmostic  coefficient  is  related  to  aw  by 

ta°-  -  -  '-w  <3) 

where  A#  is  the  molecular  weight  of  water. 

The  variation  of  G"  with  temperature  is  given 


dCr\T  -L 

9T  ~  r 


where  L  is  the  total  relative  enthalpy,  i.e.,  the  dif¬ 
ference  between  the  total  enthalpy  of  the  solution 
and  the  enthalpy  of  the  components  at  infinite  di¬ 
lution.  The  temperature  variation  of  L  is  given  by 


dT 


where  J  is  the  total  relative  heat  capacity  of  the  so¬ 
lution: 

J  =  Cpt  ~(n,Cpt+ /tj  Cp,°)  (6) 

where  Cpt  =  total  heat  capacity 

n,  =  number  of  moles  of  solvent 
n,  =  number  of  moles  of  solute 
Cpt  and  Cpt  =  partial  molal  heat  capacities  of 
water  and  solute,  respectively,  at 
infinite  dilution. 

From  total  heat  capacity  measurements  the  tem¬ 
perature  dependence  of  the  activity  and  osmotic 
coefficients  can  thus  be  found  by  integrating  eq  4 
and  5. 

The  prediction  of  activity  coefficients  of  elec¬ 
trolytes  has  recently  been  improved  by  the  specific 
interaction  equations  of  Pitzer  and  coworkers 
(Pitzer  1973,  197S,  Pitzer  and  Mayorga  1973, 
1974,  Pitzer  and  Kim  1974).  The  mean  activity  co¬ 
efficient  of  a  1:1  salt  as  a  function  of  concentra¬ 
tion  is  given  by  a  polynomial  with  second  and 
third  virial  coefficients: 

lny±  *  p+mB'1  +  miCy  (7) 

and  the  osmotic  coefficient  by 

<t>- 1  «/*+mB*  +  m,C*  (8) 

where/1,  and/*  represent  the  effects  of  long-range 
electrostatic  interactions  among  ions.  The  param¬ 
eters  B7  and  B*  are  the  second  virial  coefficients 
for  short-range  interactions  between  two  ions 
(m-m,x-x,  m-x).  The  terms  C7  and  C*  are  the 
third  virial  coefficients  for  triplet  interactions  and 
are  important  at  high  concentrations. 

By  substitution  of  the  expressions  for  y±  and 
(d- 1)  (eq  7,  8)  into  eq  2,  the  excess  free  energy  is 
now  expressed  in  terms  of  Pitzer’s  equations  by 

G"  «  nwRT(f  +  m'B^  +  (9) 

where /is  the  long-range  electrostatic  term  (P~ 


r 


2 


f*),  and  An  and  CM  are  the  short-range  interac¬ 
tion  terms  for  excess  free  energy.  The  full  forms  of 
eq  7-9  are  lengthy  and  are  displayed  in  Appendix 
A. 

The  first  and  second  temperature  derivatives  of 
eq  9  as  prescribed  by  eq  4  and  5  yield  a  temper¬ 
ature-dependent  equation  that  can  be  fit  to  total 
heat  capacity  data  to  find  the  parameters  Bm  and 
CK  at  a  given  temperature: 

Cpt  -  nwCp°  +  mCp?  + 

+  (10) 

where  and  are  the  second  derivatives  of 
the  second  and  third  virial  coefficients  and  f  is  the 
second  derivative  of  /.  The  Debye-Hiickel  slope, 

which  gives  the  relationship  of  thermody¬ 
namic  properties  to  concentration  at  high  dilution, 
is  included  in  the  /  term  (eq  7  in  Appendix  A). 
Also,  A4  depends  on  the  density  of  water  and  the 
dielectric  constant,  which  are  functions  of  temper¬ 
ature.  The  derivatives  of  A4  with  respect  to  tem¬ 
perature  have  been  calculated  by  Silvester  and  Pit- 
zer  (1978)  from  0°C  to  high  temperature.  The 
equations  for  enthalpy  and  heat  capacity  in  terms 
of  Pitzer's  equations  have  been  derived  by  Silves¬ 
ter  and  Pitzer  (1977)  and  are  given  in  Appendix  A. 


Heat  capacity  measurements 

The  instrument  used  for  measuring  specific  heat 
was  a  standard  Perkin-Elmer  Model  DSC-2.  Al¬ 
though  the  DSC-2  is  not  the  most  accurate  instru¬ 
ment  for  certain  types  of  calorimetric  measure¬ 
ments,  it  is  appropriate  for  the  present  purpose  be¬ 
cause  heat  capacities  can  be  measured  on  solutions 
well  below  their  normal  freezing  point  by  super¬ 
cooling  the  solutions.  The  temperature  that  can  be 
reached  before  freezing  varies  with  the  salt  con¬ 
centration.  Typical  thermograms  are  shown  at 
several  concentrations  of  NaCl  in  Figure  1 .  At  0.S 
tnoial  the  solution  was  cooled  to  about  250  K  be¬ 
fore  ice  formed,  as  shown  by  the  large  exotherm  at 
250  K.  As  the  concentration  increases,  first  ice  and 
then  a  eutectic  mixture  of  ice  and  solid  sodium 
chloride  precipitates  from  the  supersaturated  solu¬ 
tions  as  the  temperature  decreases.  With  increas¬ 
ing  concentration  the  ice  and  eutectic  exotherms 
move  closer  together  and  both  are  shifted  to  lower 
temperatures.  The  6.0-molal  solution  remained 
supersaturated  to  below  230  K  (-43  °Q,  which  is 
the  limit  of  temperature  control  of  the  calori¬ 
meter. 


The  DSC-2  measures  heat  flow  into  a  sample 
per  unit  time  ( dH/dt )  relative  to  a  known  standard 
(usually  sapphire).  The  specific  heat  Cp  is  related 
to  the  heat  flow  by 


dH 

dt 


mCp 


dT 
dt  ‘ 


(11) 


The  DSC-2  data  are  stored  and  processed  by  soft¬ 
ware  written  by  New  England  Digital  Corporation 
for  an  Able  60  computer  located  in  the  CRREL 
Soils  Laboratory.  The  system  collects  data  from 
the  DSC-2  and  computes  specific  heat  at  one- 
degree  intervals.  Energy  calibration  is  afforded  by 
data  from  a  sapphire  standard  obtained  over  the 


Ttmpwistur*  (K) 


Figure  1.  Thermograms  of  cooling 
scans  at  different  concentrations  of 
sodium  chloride. 


same  temperature  range.  Temperature  is  calculat¬ 
ed  by  measuring  the  freezing  temperatures  of  pure 
water  and  mercury  at  a  scan  rate  of  5°  min'1  inti  ; 
wanning  mode.  An  empty  aluminum  sample  pan 
is  run  to  establish  the  baseline.  The  samples,  emp¬ 
ty  pan  and  sapphire  were  all  run  at  a  scan  rate  of 
-10°  min'1  in  the  cooling  mode  at  the  highest  in¬ 
strument  sensity  (20  meal  s'1).  Sample  weights 
averaged  20  mg  and  were  weighed  on  a  Perkin- 
Elmer  TGS-2  thermobalance  precise  to  0.1  pg.  A 
subambient  kit  was  used  which  included  a  glove 
box  and  a  helium  purge  gas  supply.  The  system 
was  cooled  with  a  refrigeration  unit  supplied  by 
CRREL. 

The  small  sample  size  (about  0.025  cm’)  and  low 
thermal  resistance  of  the  DSC-2  sample  holder 
permits  fast  scan  rates  without  temperature  gradi¬ 
ents  and  thermal  lag  in  the  sample.  Specific  heats 
can  therefore  be  measured  to  low  temperatures  on 
supercooled  solutions.  Accuracy  is  sacrificed, 
however,  especially  at  low  concentration,  as  com¬ 
pared  with  some  adiabatic  and  flow  calorimeters, 
but  these  systems  cannot  be  used  on  samples  with 
decreasing  temperature  (adiabatic)  or  below  0°C. 
It  will  be  shown  that  accurate  high-dilution  meas¬ 
urements  are  not  needed  to  predict  adequately  the 
activities  of  saturated  solutions  below  0°C. 

Three  samples  were  prepared  at  each  concentra¬ 
tion  with  different  sample  pans.  The  sapphire 
standard  and  empty  pans  were  rerun  before  each 
set  of  three  samples.  Temperature  was  recalibrat¬ 
ed  from  thermograms  of  mercury  and  pure  water 
before  and  after  the  standards  and  samples  were 
run.  The  temperature  calibrations  were  reproduci¬ 
ble  to  about  0.05  K. 

The  heat  flow  ( dH/dt )  into  each  sample  was 
converted  to  specific  heat  by  eq  11.  The  specific 
heat  values  were  smoothed  over  temperature  by 
fitting  the  data  at  one-degree  intervals  to  second- 
order  polynomials  in  temperature.  The  standard 
error  of  a  point  from  the  smooth  curve  varied  with 
temperature  and  concentration  from  1.6ft  at  0.5 
molal  and  25°C  to  0.8ft  at  6.0  molal  and  -40°C. 

The  average  of  the  smoothed  specific  heat  and 
the  standard  error  for  the  three  replicates  at  each 
concentration  are  given  in  Appendix  B.  The  repro¬ 
ducibility  of  the  values  is  2.4ft  at  0.5  molal  and 
25°C,  and  decreases  to  0.7ft  at  6.0  molal  and 
-40°C. 

Calculations 

The  calculation  of  low-temperature  activities  of 
sodium  chloride  and  water  begins  by  determining 
the  temperature  dependence  of  the  adjustable  par¬ 


ameters  /J®,  and  C*  in  Pitzer’s  equations.  This 
dependence  can  be  found  by  fitting  activity  or 
osmotic  coefficient  data  directly  to  eq  A7,  if  such 
data  are  available.  For  most  salts,  however,  very 
few  reliable  activity  measurements  have  been 
made  below  10°C  at  high  concentration.  Silvester 
and  Pitzer  (1977)  used  apparent  heat  capacity  data 
(4Cp)  combined  with  other  available  thermody¬ 
namic  data  (enthalpy,  activity  and  osmotic  coeffi¬ 
cients)  to  extend  the  ion-interaction  equation  of 
Pitzer  from  25  °  to  300°C. 

The  heat  capacities  measured  in  this  study  can 
be  used  to  improve  the  predictions  from  25°  to 
-40°C  by  fitting  total  heat  capacity  data  to  eq  A14 
and  A15.  The  procedure  is  essentially  the  same  as 
that  described  by  Rogers  and  Pitzer  (1981)  for  so¬ 
dium  sulfate  solutions  with  the  exception  that 
total  heat  capacities  instead  of  relative  heat  capac¬ 
ities  are  used  in  this  study.  The  only  difference  is 
that  the  additional  parameter,  Cp°,  the  molal  heat 
capacity  of  pure  water,  must  be  fit  to  eq  A14  at 
each  temperature. 

First,  the  specific  heat  of  the  sodium  chloride 
solutions  is  converted  to  total  heat  capacity  of  an 
m-molal  solution  containing  one  kilogram  of 
water  by 

Cpt  «  Cp( 1000 + mM)  (12) 

where  A#  is  the  molecular  weight  of  sodium  chlor¬ 
ide.  The  values  of  Cp,°  and  Cp?  are  found  by  fit¬ 
ting  Cpt  to  eq  A14  and  treating  Cp?,  Cp°, 
and  as  adjustable  parameters.  A  linear  least- 
squares  regression  program  (Davis  1973)  was  used 
to  fit  the  data.  Below  0*C  the  Debye-Huckel  slope 
(4#,  Ah  and  A})  was  extrapolated  from  Pitzer’s 
values  above  0°C. 

The  fitted  values  of  Cp°  and  €p°  determined 
above  were  uri  to  find  the  best  values  of  the  sec¬ 
ond  temperature  derivatives  of  0®,  0°*  and  C* 
(App.  A)  at  each  temperature.  The  heat  capacity 
data  were  refit  to  eq  A14,  and  they  reproduced  the 
specific  heat  data  to  within  0.004  cal  g'1  K'\ 

The  temperature  variation  of  the  parameters 
d'&V/dT1,  P&V/dT1  and  dlC°/bTl  is  most  sim¬ 
ply  represented  by  second-degree  polynomials 
from  25  °  to  -40*C.  The  polynomial  equations  are 
integrated  by  parts  according  to  eq  A16  and  A17 
after  both  sides  are  multiplied  by  T.  The  values  of 
0®,  0°’  and  C°  and  their  first  temperature  deriva¬ 
tives  at  25  *C  were  used  as  the  integration  con¬ 
stants.  Values  of  0®,  0<l)  and  C°  at  25  #C  are  from 
White  and  Bates  (1980)  and  the  first  derivatives 
are  from  Silvester  and  Pitzer  (1977).  The  tern- 


perature  dependence  of  the  parameters  is  given  by 
the  following  equations: 


d_  AG, 
dT  T 


(17) 


0*  -  2x,  +2xi/T+x,T1/6  + 

jr4rV12+jr,r4/20  (13) 

=  2xt  +  2xi/T+x,  T*/6  + 
jr*rV12+jf,#r4/20  (14) 

C*  =  2jc,,+2x„/r+JfllrV6+ 

jr,4TV12+x„r4/20.  (15) 

The  coefficients  of  eq  13-15  are  listed  in  Table  1. 


Table  1.  Coefficients  of  eq  13-15. 

x,  =  7.66249  x,  =  4.49459  (10-’) 

x ,  *  -1.00472(10*)  x„  =  -7.86640  (10  *) 

jr,  =  -1.79001  (10  *)  x,,  a  -1.06828 

X.  -  1.15601  (10  ’)  *  1.42620(10*) 

jr,  =  -1.90346  (10  ■)  Jf„  =  2.44524  (10  *) 

x.  *  2.57468  (10*)  Jf,.  a  -1.56052(10“) 

a  -3.25771  (10*)  x„  a  2.54348  (10“) 

jr.  =  -6.58529  (10*) _ 


RESULTS 

Values  of  >±  and  aw  calculated  over  the  entire 
temperature  and  concentration  range  of  this  study 
are  listed  in  Table  2.  The  accuracy  of  these  con¬ 
stants  can  be  checked  by  independently  calculating 
ay.  from  the  thermodynamic  properties  of  ice.  In  a 
solution  of  sodium  chloride  at  the  freezing  point 

AG$"(  ice)  =  -RTlna„  (16) 

where  AG"  is  the  free  energy  change  in  going  from 
ice  to  infinite  dilution  at  some  temperature  T. 

To  calculate  the  temperature  dependence  of 
AG",  the  derivative  with  respect  to  temperature 
yields 


•°  _ 

"  “  T* 

where  A// “is  the  heat  of  solution  at  infinite  dilu¬ 
tion.  At  0®C,  A //"  =  - 1436  cal.  The  variation  of 
A //"is  found  from 

dAH? 

-if  =  ACp"  (18) 

where  A Cp°  is  the  difference  in  heat  capacity  be¬ 
tween  ice  and  pure  water: 

AQ?"(ice)  =  Cp?~Cp.  (19) 

The  value  of  Cp° is  from  the  measured  heat  capac¬ 
ities  of  water  emulsions  down  to  -38  "C  (Angell  et 
al.  1973).  The  heat  capacity  of  ice  (Cp)  is  from 
Dorsey  (1940).  The  temperature  variation  of  ACp° 
was  fit  to  the  following  equation: 

AC/»"(ice)  =  z,+ZiT+ziT2.  (20) 

Integration  of  eq  20  with  the  values  of  AH° and 
AG“for  ice  at  0"C  gives 

AGs°(ice)  =  -z.Tlnr-y  T1  - 

z,r 

— +z,T+z>.  (21> 

The  coefficients  of  eq  20  and  21  are  listed  in  Table 
3.  A  comparison  of  the  activity  of  water  calculated 
from  the  properties  of  ice  and  from  Pitzer’s  equa¬ 
tions  using  the  parameters  determined  above  is 
given  in  Table  4.  The  average  standard  deviation  is 
0.0017.  In  Figure  2  the  calculated  solubility  of  ice 
is  plotted  against  the  experimental  solubility;  the 
agreement  is  good. 

The  activities  of  components  in  a  salt  mixture  at 
low  temperature  can  be  calculated  from  the  mix¬ 
ing  equations  of  Pitzer  and  Kim  (1974)  by  using 
the  pure  salt  parameters  /3<0),  ft])  and  C°  as  deter¬ 
mined  in  this  study  and  the  mixing  parameters  6 
and  ^  for  different  ions  of  the  same  charge  tvoe. 
The  terms  0  and  are  small  at  25  °C  and  probably 
will  not  require  correction  at  low  temperature. 
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Table  3.  Coefficient* 
of  aq  20  aad  21. 


Table  2.  Valocs  of  y±  sad  a*  cakalated  below  0°C  ia  this  ctady. 


0*C 

- 

•5*C 

- 

KPC 

m 

7± 

•> 

y± 

•* 

y± 

>. 

Z,  -  3.0256(10*) 

z,  -  -1.1274(10-*) 

0.3 

0.636 

0.9039 

0.646 

0.9840 

0.633 

0.9842 

Z,  -  1.2876  (10-*) 

1.0 

0.614 

0.9679 

0.398 

0.9684 

0.578 

0.9688 

24 

0.999 

0.9347 

0.574 

0.9338 

0.346 

0.9372 

Z.  -  1.68746(100 

3.0 

0.623 

0.676 

0.8902 

0.8376 

0.392 

0.639 

0.9001 

0.8398 

0.536 

0.398 

0.9023 

0.8624 

Z,  -  -3.49782  (100 

4.0 

S.0 

0.736 

0.8122 

0.716 

0.8142 

0.670 

0.816S 

6.0 

0.S70 

0.7620 

0.829 

0.7629 

0.780 

0.7641 

_ 

■1S*C 

20*C 

_ 

23*C 

m 

7± 

I. 

y± 

*. 

y± 

*» 

Table  4.  Activity  of  water  (aw)  at  frees* 

0.3 

0.618 

0.9844 

0.600 

0.9846 

0.580* 

0.9849* 

inf  point  calculated  from  Pitzer’s  eqna- 

1.0 

0.336 

0.9694 

0.530 

0.9702 

0.501* 

0.9710* 

tioes  compared  to  values  from  AG.'Tlce) 

2.0 

3.0 

0.313 

0.316 

0.9389 

0.9049 

0.477 

0.473 

0.9408 

0.9080 

0.438 

0.426 

0.9431 

0.9117 

aad  solablHty  data. 

4.0 

3.0 

0.332 

0.619 

0.8636 

0.8193 

0.301 

0.363 

0.8693 

0.8227 

0.448 

0.304 

0.8736 

0.8266 

«w/rom 

6.0 

0.726 

0.7635 

0.667 

0.7673 

0.603 

0.7693 

T 

Pitzer’s  Solubility 

rc> 

equation  data 

Difference 

-30*C 

- 

33*C 

- 

40 *C 

-3.48 

0.9666  0.9668 

0.0002 

Iff 

y± 

r± 

•* 

7± 

•» 

-5.17 

0.9315  0.9510 

-0.0003 

-6.32 

0.9414  0.9403 

-0.0009 

0.3 

0.337* 

0.9832* 

0.513* 

0.9836* 

0.303* 

0.9860* 

-8.52 

0.9222  0.9205 

-0.0017 

1.0 

0.469* 

0.9720* 

0.434* 

0.9731* 

0.397* 

0.9744* 

-9.41 

0.9130  0.9126 

-0.0024 

2.0 

0.396* 

0.9438* 

0.332* 

0.9489* 

0.307* 

0.9326* 

-11.04 

0.9007  0.8982 

-0.0023 

3.0 

0.377 

0.9160 

0.327* 

0.9210* 

0.278* 

0.9268* 

-14.77 

0.8682  0.8660 

-0.0022 

4.0 

0.393 

0.8788 

0.317 

0.8848 

0.282* 

0.8918* 

-16.21 

0.8362  0.8539 

-0.0023 

3.0 

0.443 

0.8313 

0.380 

0.8368 

0.319 

0.8432* 

-18.73 

0.8346  0.8331 

-0.0013 

6.0 

0.336 

0.7718 

0.467 

0.7748 

0.397 

0.7784 

-20.83 

0.8129  0.8160 

0.0031 

•  Values  ffltfwl  temperature  range  of  heat  capacity  equations. 


NeCi  Coacsatrtioa  (motomy) 

Figure  2.  freezing  point  of  ice  in  sodium  chloride  solutions.  The  cunt 
is  the  experimental  solubility  pom  Links  (1958).  The  points  are  pom  this 
study. 
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CONCLUSION 

This  work  demonstrates  that  total  heat  capacity 
data  can  be  used  to  improve  calculations  of  activi¬ 
ties  of  salts  at  freezing  temperatures.  Now  that  the 
usefulness  of  the  method  has  been  shown,  addi¬ 
tional  salts  will  be  investigated  so  that  mineral 
equilibria  can  be  determined  in  common  brine 
mixtures.  In  mixtures  with  many  possible  phases, 
the  phase  and  solution  composition  can  be  found 
from  the  temperature  dependence  of  the  free 
energy  of  all  species  that  can  be  present  by  using 
one  of  the  free  energy  minimization  methods  that 
are  in  use  (Storey  and  Van  Zeggeren  1964,  Harvie 
and  Weare  1980). 
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APPENDIX  A:  PITZER’S  EQUATIONS  FOR  THE  ACTIVITY, 
EXCESS  FREE  ENERGY,  ENTHALPY  AND  HEAT  CAPACITY 

OF  A  SINGLE  ELECTROLYTE 

r 

Osmotic  coefficient 

i  I  r  -  1 4  /,/2  ,  m2vmvto  n*  .  (»’»*’x)3/2  _ 

0  1  |ZmZxl'4*1  +  1<2/»/2  »  &mx  +  m  2  , 

(Al) 

r 

B*^  =  0® + /8(,)  exp(  -  a/,/2) 

(A2) 

*■ 

Activity  coefficient 

r 

+  m2v*  /%,+„■ 

j;  mx  p  mx 

(A3) 

Wnx  -  20<o,+  ^  [l-(I  +  a/,/J-  y7  exiK-a/'^J 

(A4) 

r 

<3*  =  3/2 

(AS) 

Excess  Gibbs  free  energy 

AG**  =  rtyfVmR  T(1  —  0+ln7±) 

(A6) 

r 

AG“  =  /iwR7j-A*  ^bKl  +  1.2fl^  +  m>2rarsR^-t-m*2(rBr]l)>/2Qn] 

(A7) 

Bmx  -  0<O)  +  ^  U  -0  +«/,/J)  exp(-«/,/J)] 

(A8) 

- 

■ 

c 

t-mx  -  2 

(A9) 

Relative  enthalpy 

**■* 

-v. 

L  -  mv|ZmZx|^  In(l  + 1 .2Il/1) - 2vm»xRrim'Bix+2(rmrJi/2RT'm>C‘iiX 

(A10) 

_ 

BL  «  (JV^M  -(^  +  (^5l(J/)II‘(I+tt/,,J) «P<-a/,/J)l 

(All) 

CL  3  /3C«w\  1 

L,M  \ar/p  2\dr/p 

(AI2) 

(A13) 

- 

9 

•* 

f 

y 

> 

Relative  heat  capacity 


y  =  mHZmZ*l^  \n(l  +  l.2Il^-2>m*xRrrn1B{m-KpmvJi*Rrm'CJmx  (A14) 


J  =  Cpt-nwCP°-mCpf  (A15) 


“  Var*  /p.i  +  T\dTk.  i 

(A16) 

lr0«»\  2  Id#®]  [ 

”  \ar  Jr  +  T\dTj p  +| 

Va^»\  .  2  1_2_ 

.yar* /p  r  \ar/rfP/ 

x  (1  -(1  +ali/3)  exp(-a/,/2)] 

w  _  .  2  ldC™\ 

mx  \ar  / p  T\dT  ) p 

(A17) 

ft. 

SET 

«S_Pt 

ii 

(A18) 

Nomenclature 

Debye-Hiickcl  slope  for  the  activity  and  osmotic  coefficient 
Ah  Debye-Huckel  slope  for  enthalpy 
Aj  Debye-Huckel  slope  for  heat  capacity 
Zm,  zx  charge  on  cation  and  anion,  respectively 

I'm,  vx  number  of  cations  and  anions,  respectively,  per  formula  weight  in  a 
strong  electrolyte  mx 
a  2.0,  a  constant 

/  ionic  strength,  I  =  xAm  vxzi) 

nw  mass  of  water  (kg) 
m  molality  (mol  salt- kg  H,0'*) 

Cpt  total  heat  capacity  (cal  K'1  mol'1) 

Cp,°  molal  heat  capacity  of  pure  water  (cal  K'1  mol'1) 

Cp°  partial  molal  heat  capacity  of  solute  at  infinite  dilution  (cal  K~‘  mot'1) 
ff®,  ft"  and  C£,  fitted  parameters 
AG"  excess  free  energy 
R  gas  constant 
T  temperature  (K) 

L  total  relative  enthalpy 
J  total  relative  heat  capacity 


APPENDIX  B:  AVERAGE  VALUES  OF  SPECIFIC  HEAT  AND  STANDARD  ERRORS 
FOR  THREE  RUNS  AT  EACH  CONCENTRATION 


CofomtnUton  23*C  20*C  J3*C  10*C  5*C  0*C  -S'C 

(mol  kg-')  Cp  EE  Cp  EE  Cp  S.E.  Cp  EE  Cp  EE  Cp  EE  Cp  EE 


0.502 

1.000 

0.026 

1.001 

0.026 

1.002 

0.0264 

1.003 

0.0266 

1.003 

0.0267 

1.0038 

0.0264 

1.005 

0.026 

1.000 

0.9400 

0.0039 

0.9490 

0.0030 

0.9504 

0.0036 

0.9506 

0.0034 

0.9504 

0.0031 

0.9499 

0.0029 

0.9490 

0.0027 

1.995 

0.9076 

0.0106 

0.9001 

0.011 

0.9001 

0.0111 

0.9075 

0.0110 

0.9064 

0.0106 

0.9048 

0.0090 

0.9024 

0.009 

2.550 

0.0015 

0.0040 

0.8030 

0.003 

0.8838 

0.0022 

0.8837 

0.0017 

0.8828 

0.0013 

0.0812 

0.0013 

0.8770 

0.0042 

3.075 

0.0549 

0.0175 

0.0557 

0.019 

0.8558 

0.0193 

0.8555 

0.0196 

0.8546 

0.0196 

0.053 I 

0.0193 

0.8510 

0.01S7 

4.096 

0.0310 

0.0050 

0.8323 

0.00S9 

0.8329 

0.0066 

0.8328 

0.0072 

0.8319 

0.0077 

0.8303 

0.0078 

0.8280 

0.0079 

5.027 

0.0000 

0.0147 

0.8111 

0.0150 

0.8126 

0.0151 

0.0133 

0.0150 

0.0133 

0.0147 

0.0125 

0.0140 

0.8110 

0.0132 

6.002 

0.7073 

0.0120 

0.7891 

0.0115 

0.7903 

0.0109 

0.7909 

0.0104 

0.7909 

0.0099 

0.7903 

0.0094 

0.7891 

0.0008 

Conctntrwrkm 

-10*C 

-15  *C 

-30*C 

-35*C 

-iO'C 

-35  "C 

-40  *C 

(mol  kg-') 

EE- 

Cp  EE- 

Cp  S.E. 

Cp  S.E. 

Cp  S.E. 

<5  EE- 

Cp  EE- 

0.502 

1.0051 

0.0258 

1.0059 

0.0258 

1.0068 

0.0258 

1.0073 

0.025 

1.0079 

0.0250 

1.0006 

0.0246 

1.0092 

0.0244 

1.000 

0.9477 

0.0026 

0.9460 

0.0025 

0.9440 

0.0027 

0.9416 

0.0030 

0.9308 

0.0055 

0.9357 

0.0041 

0.9322 

0.0050 

1.995 

0.8997 

0.0077 

0.8963 

0.0063 

0.8925 

0.0050 

0.0000 

0.0042 

0.8831 

0.0048 

0.8776 

0.0068 

0.8715 

0.0096 

2.558 

0.8754 

0.0018 

0.8713 

0.0022 

0.8665 

0.0028 

0.8600 

0.0035 

0.8636 

0.0119 

0.8470 

0.0055 

0.8389 

0.0069 

3.075 

0.8484 

0.0170 

0.0452 

0.0166 

0.0415 

0.0151 

0.0372 

0.0135 

0.8323 

0.0117 

0.8269 

0.0099 

0.8209 

0.0006 

4.096 

0.8250 

0.0070 

0.8213 

0.0075 

0.0168 

0.0072 

0.8117 

0.0067 

0.8058 

0.0062 

0.7992 

0.0058 

0.7918 

0.00S7 

5.027 

0.0007 

0.0120 

0.8056 

0.0106 

<r.0018 

0.0009 

0.7973 

0.0071 

0.7919 

0.0054 

0.7859 

0.0044 

0.7790 

0.0053 

6.002 

0.7873 

0.0002 

0.7050 

0.0076 

0.7020 

0.0071 

0.7784 

0.0066 

0.7743 

0.0061 

0.7696 

0.0056 

0.7642 

0.0054 

A  facsimile  catalog  card  in  Library  of  Congress  MARC  format  is  repro¬ 
duced  below. 
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